Hydrogen Release from Sodium Alanate 
Observed by Time-resolved Neutron Backscattering 
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Innermolecular motion in Na3AlH6 gives rise to a Lorentzian spectrum with a wavenumber- 
independent width of about 1 /ieV at 180° C, which is probably due to rotation of AIHe tetrahedra. 
There is no such quasielastic line in NaAlIrU or NaH. Based on this finding, time-resolved measure- 
ments on the neutron backscattering spectrometer SPHERES were used to monitor the decompo- 
sition kinetics of sodium alanate, NaAlEU — > NaaAIHe — > NaH. Both reaction steps were found to 
be accelerated by autocatalysis, most likely at the surfaces of NasAlHg and NaH crystallites. 
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I. INTRODUCTION 

Sodium alanate NaAlH4 is widely studied as a model 
system for hydrogen storage PHI]- While hydrogen ex- 
change in pure bulk alanate is not reversible, hydrogen 
desorption and absorption kinetics has been improved de- 
cisively by doping 0, @ and nanostructuring [6-8] . Fur- 
ther improvement of the hydriding/dehydriding kinetics 
requires a detailed understanding of the reaction mecha- 
nism. 

In this work, we introduce time-resolved neutron 
backscattering as a new tool to monitor the solid-state 
reactions that are responsible for hydrogen release. Meas- 
urements are carried out on pure NaAlH4 and its decom- 
position products Na3AlH 6 and NaH; one explorative ex- 
periment has been undertaken on material doped with 
TiCl3, cycled under hydrogen, and quenched in the hy- 
drogenated state. 



A. The Hydriding/Dehydriding Reactions 

Hydrogen storage in sodium alanate is based on the 
two reactions 



NaAlH 4 O ±Na 3 AlH 6 + |A1 + H 2 



Na 3 AlH 6 o 3NaH + Al+§H 2 . 



(1) 
(2) 



A third reaction, the decomposition of NaH, has no prac- 
tical importance, because it only occurs above 425°C |9j. 
Therefore, a theoretically reversible capacity of 5.6 wt.% 
hydrogen is available. The equilibrium temperatures of 
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([T]) and J2]) depend only little on doping, but much on 
pressure [10]. At 1 bar, they are about 30°C and 110°C, 
respectively, with uncertainties of several degrees [ll], 
Fig. 2]. A thermodynamic argument suggests that at 
particle sizes below 50 nm, NaAlH/t may decompose in a 
single step into NaH, Al, and H 2 [13]. On the other hand, 
time-resolved neutron diffraction during hydrogen release 
from NaAlH 4 with a particle size of 110 nm clearly shows 
a transitory state rich in NaaAlHg [13 ]. 

In undoped bulk NaAlH4, reaction rates for hydrogen 
release are very low, unless the temperature approaches 
the melting point for which values between 175°C and 
190°C are reported tlMl]. In pure ball-milled NaAlH 4 
at 170°, it takes at least 8 h to complete the desorption 
step ([lj. When doped with Ti clusters, for example, this 
step takes only 300 s at 150°C [H. 

From studies on doped material, it can be concluded 
that the reactions (p}, ([2]) take place at the boundary be- 
tween Al particles and NaH, NaaAlHg, or NaAlH 4 phases 
[ItJ and that the rates of these reactions are limited by a 
nucleation/growth process [l~8| or/and by the transport 
of heavy, Al- or /and Na-based species [3, EH. This is sup- 
ported by H/D exchange experiments showing that the 
rate-limiting step of the hydriding/dehydriding process 
is neither H 2 dissociation nor H diffusion [2(J [2l[ . 

Quite different mechanisms have been proposed to ex- 
plain the effect of titanium doping: (i) Ti is a surface 
catalyst 0, [221423 ]; (ii) Ti creates Na + vacancies in 
the bulk structure, thereby facilitating hydrogen diffu- 
sion Jill [25|, [26[ or/and the transport of heavier species 
Il9j ; (iii) Ti weakens the Al-H bond, promoting the 
removal of H 2 [27j ; (iv) TiCl3 is a grain refiner preventing 
growth of Al and NaH particles [l~3| ; (v) Ti is an initiator 
of Al nucleation [|J] . 
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B. Quasielastic Neutron Scattering Studies 

We are aware of two previous quasielastic neutron scat- 
tering measurements. One measurement using a triple- 
axis spectrometer on NaAlH4 at 150°C indicates that 
a small fraction (< 10%) of H atoms participates in a 
fast localized process (almost q independent Lorentzian 
halfwidth r ~ 70 /ieV, corresponding to a characteris- 
tic time of r = h/T ~ 10 ps) [28[. Doped and undoped 
NaAlH4 and NaaAlHg were investigated using the Jiilich 
backscattering spectrometer BSS [29|, [3(J. In doped 
NaAlH 4 at 117°C, about 0.5% of the hydrogen were 
found to be mobile, too little to make reliable statements 
about its motion [2{|. In doped Na 3 AlH 6 at 77° C and 
in undoped NaAlH 4 at 117°C, some quasielastic broad- 
ening was observed although the Lorentzian width T was 
considerably smaller than the resolution fwhm of 0.8 fieV 
30] . As no spectral fits were shown, it is not possible to 
assess the conclusion of broadening being due to jump 
diffusion. 

In the meantime, Forschungszentrum Jiilich has re- 
placed the BSS by the new backscattering spectrometer 
SPHERES 31]. Compared to BSS, SPHERES has much 
higher count rates, better resolution (0.65 /xeV), a wider 
dynamic range, and a better signal-to-noise ratio. It is 
now possible to measure a meaningful spectrum within 
a fraction of an hour. This enables us to observe micro- 
scopic dynamics during hydrogen release in real time. 

In this work, we will use neutron backscattering to 
monitor the reactions (fl) and @. Using the quasielastic 
amplitude as a proxy for the amount of NaaAlHg in the 
sample, we will show how this compound builds up in 
the first reaction step and disappears in the slower sec- 
ond step. Our results indicate that both reaction steps, 
after starting very slowly, are accelerated significantly 
when autocatalysis on the surface of NasAlHg and NaH 
crystallites sets in. Preliminary data for TiCbj-doped 
NaAlH4 suggest that the dopant acts like a nucleation 
center, accelerating the initial, non-autocatalytic phases 
of both reactions. 



II. SAMPLE PREPARATION 

In this paper, we report on four sample materials: 
NaH, Na 3 AlH 6 , NaAlH 4 , and NaAlH 4 (a8a). All sam- 
ples were prepared in an argon-filled glove box equipped 
with a recirculation system to keep the water and oxygen 
concentrations below 1 ppm. The raw chemicals NaH 
(95%, Sigma Aldrich), NaAlH 4 (96%, Albemarle), and 
TiCLj (99.999%, Sigma Aldrich) were used as received. 
Mechanical milling was carried out in a Fritsch P6 plan- 
etary mixer/mill at a rotation speed of 600 rpm using a 
silicon nitride vial and balls with a ball to powder weight 
ratio of about 20:1. The vial was filled and sealed in the 
glove box under argon atmosphere. 

The samples were prepared as follows: 

NaH: Powder as received. 



NaaAlHg: Obtained by mechanical alloying of NaH 
and NaAlH 4 at a molar ratio of 2 to 1. The powder was 
milled for 20 h in an argon atmosphere as described in 
Refs. @, 0- The product was analyzed by X-ray diffrac- 
tion. 

Pure NaAlH 4 : Powder as received, ball-milled for 30 
minutes. 

NaAlH 4 (a8a): Doped with TiCl 3 , cycled eight times 
under hydrogen, and quenched in the hydrogenated state. 
This sample was obtained by ball milling 2 g of NaAlH 4 
and 285 mg of TiCL;, resulting in 5 mol.% Ti doping on 
the basis of TiCl3. In order to avoid any increase in the 
temperature, the 3 hours of milling were divided into 30 
minutes of milling and 10 minutes pause with 5 repeti- 
tions. Directly after milling, the sample was first dehy- 
drogenated, then cycled eight times under hydrogen with 
hydrogen absorption (100 °C, 100 bar) and desorption 
(150 °C, 0.3 bar) conditions, and finally quenched in the 
hydrogenated state. Absorption and desorption of hydro- 
gen were carried out in a carefully calibrated modified 
Sieverts apparatus. A more detailed description of the 
apparatus, the reactor, and the absorption/ desorption 
procedure can be found elsewhere [la, |32| • 

Unfortunately, lack of beamtime prevented us from 
measuring another doped sample, freshly prepared with- 
out hydrogen cycling. 

The samples were filled into flat, top- loading 30 x 40 x 
0.5 mm 3 Al cells and sealed with Al wire in the glove 
box. In the course of this study, we were concerned about 
the sample material, while releasing hydrogen, possibly 
reacting with the cell. Therefore, we performed an ad- 
ditional measurement with NaAlH 4 powder filled into a 
pocket made of Ag foil to isolate the sample from the Al 
walls. Direct comparison of the neutron spectra revealed 
no difference. 



III. MEASUREMENTS AND DATA ANALYSIS 

Since the samples are rich in hydrogen, the neutron 
scattering cross section is dominated by incoherent scat- 
tering by hydrogen. 

Measurements were performed using the neutron 
backscattering spectrometer SPHERES of the Jiilich 
Centre for Neutron Science (JCNS) at the neu- 
tron source FRM II (Forschungs-Neutroncnqucllc Heinz 
Maier-Leibnitz) in Garching, Germany [3l|. To obtain 
reasonable statistics within short time slices we did not 
use the instrument's full energy range of ±31 /xeV; in- 
stead, the Doppler velocity amplitude was set to 1.3 m/s, 
giving access to a window of ±8.6 fieV. 

The sample was mounted in a cryofurnace and oriented 
at 135°. In this standard geometry, self-shielding within 
the cell limits the range of useable scattering angles to 
26 < 125°. This results in nine large-angle detectors with 
wavenumbers q = 0.6 . . . 1.8 A -1 . The small-angle detec- 
tors at q = 0.25... 0.46 A" 1 are not in exact backscat- 
tering and therefore have comparatively bad resolution. 
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The raw data reduction was carried out with SLAW 
331. and further analysis was performed with FRIDA 
34} . Spectra were normalized to the integral intensity at 
about room temperature. The room temperature meas- 
urements also served as resolution function. 

Deviating from common practice, we did not subtract 
an empty-cell measurement. Instead, cell scattering was 
taken into account in the fitting procedure. In the energy 
range of SPHERES, scattering from Al just consists in a 
weak elastic line. Besides, there is a flat background due 
to various imperfections of the instrument. This back- 
ground was determined from the baseline of the resolu- 
tion spectrum, and included as a fixed component in all 
fit models. In fitting, squared deviations were weighed 
with the reciprocal standard deviation, and theoretical 
expressions are convoluted with the measured resolution 
function (minus the aforementioned baseline). 

The data were analyzed in an heuristic and itera- 
tive way, trying different procedures and testing differ- 
ent models before a consistent overall picture emerged. 
In the following sections, data analysis will be decou- 
pled from interpretation. In Sect. lIVl spectral line shapes 
are analyzed, and simple fit functions are introduced. In 
Sect. El fits are used to extract information from time- 
resolved measurements. Finally, in Sect. I VII the obser- 
vations are interpreted physically. 



IV. SPECTRA 

The simplest fit function that works for the entire data 
set is 

S(q,u;) = b + asS(Lj) + ^2ai(q)C(u} ) T i ). (3) 



have been rescaled to account for the temperature de- 
pendence of the elastic scattering (the Debye- Waller fac- 
tor, or more precisely the Lamb-Mossbauer factor, since 
we are talking about self correlations measured by inco- 
herent scattering). As a result, the spectra coincide per- 
fectly, demonstrating the absence of quasielastic broad- 
ening. 
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FIG. 1: Neutron scattering spectra of NaH measured using 
the backscattering spectrometer SPHERES. Here and in the 
following figures, spectral data are averaged over the q range 
0.9... 1.8 A' 1 . Only in this fi gure, spectra are normalized to 
unit area to account for the temperature dependence of elastic 
scattering. There is no temperature-dependent quasielastic 
scattering; the spectra represent just the instrumental resolu- 
tion function. 



It consists of a flat background, an elastic S line, and 
zero, one, or two Lorentzians 



r 



t 



7T U 2 + r 2 ' 



(4) 



The flat background is kept constant at the value de- 
termined at room temperature, except when the overall 
scattering decreases with time because of hydrogen de- 
sorption. 

Unconstrained fits indicate that the Lorentzian 
linewidths do not vary systematically with q, whereas 
the amplitudes do. The amplitudes of the Lorentzians 
increase with q before they reach a plateau at about 
1 A -1 . Therefore, in our final analysis of spectral line- 
shapes and time series, we average S(q, ui) over the q 
range 0.9. . . 1.8 A -1 to obtain a strong quasielastic sig- 
nal with good statistics. 



A. NaH 

Spectra of the final reaction product NaH were meas- 
ured at 38, 117, and 177°C. In Fig.[TJ these three spectra 



B. Na 3 AlH 6 

The richest set of spectra was obtained for the inter- 
mediate reaction product NasAIHg. Spectra of a freshly 
prepared sample were measured up to 177°C, where de- 
composition set in. Kinetic aspects will be discussed be- 
low (Sect. EJ. 

Two Lorentzians are needed to describe the quasielas- 
tic scattering. To reduce the number of free parameters 
and to avoid unwanted degeneracies (discussed recently 
in another SPHERES data analysis [35|]), amplitudes are 
fixed at temperature-independent values a\ — 0.5 and 
a-2 = 0.07. Excellent fits are obtained, as shown in Fig. [5] 

To investigate the q dependence of the amplitudes 
&i,2(q), we impose the linewidths Pi,2(T) as determined 
from the fits to the q averaged spectra. Results are shown 
in Fig. [3] As anticipated, the 011,2(3) are temperature- 
independent within experimental accuracy. Their q de- 
pendence is fitted reasonably well by the simple expres- 
sion 



q 



'1.2 ~~2 



+ K 2 



(5) 
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FIG. 2: Neutron scattering spectra of Na 3 AlH 6 . The 27° C 
data represent the resolution function. The other data are fit- 
ted according to Eq. <(3j with a delta line and two Lorentzians 
(here and in the following figures, all spectral fits are under- 
stood to be convoluted with the measured resolution). The 
dashed (dotted) line shows the Lorentzian 1 (2) for 117°C 
separately and in unconvoluted form. 




q (A 1 ) 

FIG. 3: Wavenumber dependence of the amplitudes ai,2 of the 
two Lorentzian components found when fitting the NaaAIHe 
spectra with Eq. @. Solid lines: Fits with Eq. 



that interpolates between a q 2 dependence for q <C k and 
a constant asymptote for g>K. 

Fig. 2] shows that the temperature dependence of the 
linewidths Ti^{T) is compatible with an Arrhenius law 
Ti = exp(— Ai/T). The activation energies are A% = 
38 kJ/mol and A2 — 23 kJ/mol; the prefactors are 
indicated in the figure. In estimating these parame- 
ters, some data points were excluded: linewidths below 
0.1 ^eV were dropped, because they are too far below 
the instrumental resolution; data at 177°C had to be ex- 
cluded, because the scattering signal was not stationary, 
so that the imposition of a fixed Lorentzian amplitude 
would not have been justified. 




□ 



2.2 2.4 2.6 2.8 
1000 K/T 

FIG. 4: Linewidths Ti^T) of the two Lorentzians used to 
fit the spectra of NaaAlHe. Only data shown as full symbols 
have been used to determine the Arrhenius laws shown as 
straight lines. The dashed line is obtained when proton NMR 
data HH are taken into account (see Fig. Illl r? = 2.3 meV, 
AS = 3600 K). 



C. NaAlH 4 
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FIG. 5: Neutron scattering spectra of NaAltL for various 
time intervals, after heating to 177° C. 

Time-resolved measurements were performed with four 
fresh NaAlH4 samples. After measuring the resolution at 
room temperature, the sample was quickly heated to 170, 
180, or 185°C, and kept there at constant temperature 
for 2 days or more. Initially, no quasielastic scattering 
is observed. It takes several hours before a quasielastic 
signal appears above the wings of the resolution function. 
The signal continues to grow, until a maximum is reached 
about 30 h after the start of the measurement. Then, the 
quasielastic intensity decreases and about 65 h after the 
start of the measurement we are left with purely elastic 
scattering, which however is only half as strong as in the 
beginning. 

To prepare for the detailed investigation of this kinet- 
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ics in Sect. [V] we need to characterise the quasielastic 
lineshapes. To do so, the data are averaged over time 
slots of 5 h (Fig. [5]) . The minimal fit function that de- 
scribes the entire data set consists of a flat background, 
a S line, and one Lorentzian. In contrast to the above 
analysis of NaaAlHg, the background must not be kept 
constant, because the total scattering intensity decreases 
during the experiment. 

The Lorentzian linewidth at 180°C and for large q is 
about 1 /ieV. This is perfectly compatible with the Ar- 
rhenius law number 1 of Fig. @] The slightly smaller 
and larger linewidths at 170 and 185°C are also com- 
patible with this Arrhenius law. We therefore attribute 
the quasielastic scattering observed during the decom- 
position of NaAlH4 entirely to the Lorentzian 1 that is 
dominant in NaaAIHg. In contrast, in the entire NaAlH4 
time series we find no trace of Lorentzian 2. 



D. NaAlH 4 (a8a) 
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FIG. 6: Neutron scattering spectra of the TiCl3-doped 
NaAlHU (a8a) sample. A reduced w scale was chosen to em- 
phasize the weak quasielastic scattering. 

Due to a lack of beamtime, only one time- resolved mea- 
surement was performed with the TiCl3-doped NaAlH 4 
(a8a) sample. After measuring the resolution at room 
temperature, the sample was heated in steps up to 118°C, 
where the onset of hydrogen release was noticed. During 
20 h, the scattering intensity decayed by more than 50%. 
Before loosing the hydrogen completely, the temperature 
was increased to 131, 137, 167, 182°C for rather short 
measurements. 

Fig.[5]shows the weak quasielastic scattering at 118°C. 
The fit consists of a free background, a S line, and a 
Lorentzian, with a fixed width of 0.20 /xeV derived from 
the Arrhenius description (Fig. 0|) of Lorentzian 1 of the 
NasAlHg analysis. 



V. KINETICS 

The results of the spectral analysis shall now be used 
for a quantitative description of the time-resolved meas- 
urements. Originally, one spectrum was saved every 5 
minutes. In our analysis, we binned them into blocks 
of 20 min. While 20 min spectra are quite noisy, it is 
perfectly possible to perform a full spectral fit. The well- 
measured resolution function is taken into account just as 
in the last section. This fitting which can be completely 
automatized must, of course, not be construed as a model 
validation. But given a valid model, it is an efficient way 
of parameter extraction. The only condition is that fit 
parameters must not be close to degeneration. There- 
fore, it is of paramount importance that fixed Lorentzian 
linewidths as obtained from the Arrhenius laws of Fig. [U 
are imposed. 

A. Na 3 AlH 6 




t (h) 

FIG. 7: Time-resolved backscattering measurement of 
Na3AlH6. Left scale: Experimental temperature sequence. 
Right scale: Total backscattering intensity, elastic intensity, 
and Lorentzian amplitudes. 

Our time-resolved measurement of NasAlHg is sum- 
marized in Fig. [7] The applied temperature sequence is 
shown along with the total backscattering intensity I (the 
integral of S(q,u) over the experimental energy range 
—8. 6... 8. 6 fJ.eV) and the outcome of the fits, namely, 
the elastic intensity ag and the two Lorentzian ampli- 
tudes 0.1,2- 

When heating step by step from room temperature 
to 152°C, I and as exhibit a parallel decrease: This 
is, at least qualitatively, the expected evolution of a 
Lamb-Mossbauer factor. Reliable amplitudes a± t 2 can 
only be determined from 97°C onwards. In the range 
97 . . . 152°C, both Lorentzian amplitudes are basically 
constant, which validates our fitting method with fixed 
linewidths Ti t %. 

Shortly before heating to 147°C, there is a transient, 
complementary excursion in as and a±, which we cannot 
explain. Another strange feature is the complementary 
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step in as and 02 after heating to 177°C. It is related to 
the outlier r 2 (177°C) in Fig.|H indicating a limitation of 
our fit. 

We measured about 6 h at 177°C. During this time, the 
quasielastic amplitude a\ decreased by 0.14. This was 
partly compensated by an increase in elastic scattering 
of 0.05 so that the total backscattering decreased only by 
0.06. 

At this point, we decided not to wait for complete de- 
composition but to cool back to 117° C for comparison 
with spectra measured before the excursion to higher 
temperatures. Comparison revealed a decrease Aai — 
—0.15, Ad2 = —0.03, partly compensated by an increase 
Aas = 0.10, resulting in a decrease of the total backscat- 
tering of AI = -0.09. 

B. NaAlH 4 




25 50 75 100 
t (h) 

FIG. 9: Quasielastic amplitude a\ in four experiments on 
NaAlHU versus time. Time scales are arbitrarily shifted to 
keep the curves separated. "Ag" designates a sample that 
was isolated from its Al container by a silver foil. Solid lines 
are fits with the same kinetic model as in Fig. [8] 
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FIG. 8: Transformation of NaAlH 4 at 185°C: Total backscat- 
tering intensity I and fitted amplitudes as and ai [Eq. ©] 
versus time. Solid lines are fits with a kinetic model described 
in Sect. EUl 

As already shown in Fig. [51 quasielastic scattering 
slowly emerged after heating NaAlH4 to a temperature 
in the range of 170. . . 185°C, only to disappear some 10 h 
later. For a closer analysis, all 20 min spectra were fitted 
with a free background, a 5 line, and one Lorenztian of 
fixed width Ti(T) taken from the Arrhenius fit of Fig. 2) 

Fig. [8] depicts the total backscattering intensity / and 
the amplitudes as, a\ as function of time, along with a fit 
that will be discussed below in Sect. lVI Bl The quasielas- 
tic amplitude a\ reaches a maximum after about 18 h, 
before it decays slowly. In the long term, this decay is 
exponential, with a half life of about 10 h. The elastic 
intensity as decreases to less than 1/6 of its initial value, 
before recovering up to about 1/4 of the room tempera- 
ture value. Total backscattering decreases continuously. 

This time-resolved measurement was carried out four 
times at three different sample temperatures of 170, 180, 
and 185°C. Fig. [5] shows the time dependence of a\ for 



all four runs. With increasing temperature, the peaks 
are somewhat sharper; in particular, the initial slope is 
much steeper. There is some fluctuation in the maximum 
value of ai, but altogether the curves demonstrate a very 
statisfactory reproducibility of measurements and data 
analysis. 



C. NaAlH 4 (a8a) 
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FIG. 10: Time-resolved backscattering measurement of 
TiCl3-doped NaAlH4 (a8a). Left scale: Experimental temper- 
ature sequence. Right scale: Total backscattering intensity, 
elastic intensity, and Lorentzian amplitudes. Solid lines are 
fits with the same kinetic model as in Fig. [8] to be discussed 
in Sect. lyTCl 

For the TiCbj-doped NaAlH4 (a8a) sample, we have 
so far undertaken one measurement only, with an im- 
provized, unsystematic temperature sequence, as shown 
in Fig. [10] Quasielastic scattering seems to set in be- 
low 80°. However, fits at such low temperatures are not 
reliable, because the linewidth Ti(T) is far below the 
instrumental resolution. Once the sample is heated to 
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118°C, the quasielastic amplitude reaches its maximum 
within less than 2 h. Then, a slow decay sets in, with 
a half life of about 14 h. The temperature excursion to 
167 and 182°C is accompanied by a dip in a\: It is not 
clear whether this is physical; it depends critically on the 
Arrhenius fit used for imposing Ti(T). 



VI. INTERPRETATION 
A. Quasielastic Scattering by Na3AlH6 

The Lorentzians used to fit the quasielastic scatter- 
ing have temperature-independent amplitudes a(q) and 
wavenumber-independent linewidths T(T) in a good first 
approximation. This is the signature of localized pro- 
cesses. We see no scattering by long-ranged diffusion: 
Most probably, diffusion of H2 is too fast and diffusion 
of heavier species is too slow to be observed within the 
dynamic window of SPHERES. 

Lorentzian 1 is due to internal motion of NaaAlHg: 
It is dominant in the quasielastic scattering of a freshly 
prepared sample, and found as a transient during the 
decomposition of NaAlH4, as expected from the two-step 
reaction formula (CO), ©. 

In contrast, Lorentzian 2 only is a small contribution to 
scattering of the fresh NasAlHe sample; it is not visible 
during decomposition of NaAlH4. Its origin must be left 
unresolved. 

The observed q dependence of the quasielastic ampli- 
tudes 01,2) approximately described by Eq. (pjj, is not 
compatible with elementary jump models that would re- 
quire an oscillatory q dependence a oc [1 — io(dq)] (jump 
length 0(d), spherical Bessel function j ) [37fl . One pos- 
sible explanation could be that each Lorentzian results 
from more than one internal mode, with similar frequen- 
cies, but different jump lengths so that ([5]) results as an 
average of 1 — jo (dq) over several d. 

Rotational motion in NaaAlHg below room tempera- 
ture was studied previously by NMR [36j, |38(. Results 
indicated thermally activated rotational jumps of AlHg 
octahedra around C4 axes. Jump frequencies, described 
by Arrhenius laws, are indicated in Fig. 1111 The slow 
C4XY rotations of Ref. [36| were found to admit activa- 
tion energies between 20 and 30 kJ/mol. Choosing the 
latter value, the Arrhenius law (the dashed line in Figs. 
2] and [TTj) extrapolates quite well towards the linewidths 
of our Lorentzian 1. This further supports our interpre- 
tation of the Lorentzian 1 being due to internal rotations 
in Na 3 AlH 6 . 



B. Decomposition of NaAlH4 

The increase and decrease of quasielastic scattering 
(Figs. 03 [5]) reflect the production and consumption of 
NasAlHg in the two-step reaction (fTJ, ([2]), respectively. 
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FIG. 11: Rotation frequencies in NasAlHg as determined by 
quasielastic neutron scattering (QENS) and by nuclear mag- 
netic responance (NMR). The QENS data are the same as in 
Fig. 01 The data point NMR 1 and the corresponding Arrhe- 
nius law was taken from [3S|. The Arrhenius law NMR 2 is 
from [36( |. The dashed line is the same as in Fig. [4] it suggests 
that NMR 2 and QENS 1 are one and the same process. 



For better readability, we abbreviate these reaction equa- 
tions as A — >• B — >• C, where A stands for NaAlH4, B for 
1/3 Na3AlH 6 , and C for NaH. Side products of the re- 
actions are not denoted, because they are irrelevant for 
the observed spectra: Al has a negligible neutron cross 
section, and H2 leaves the sample so rapidly that it does 
not contribute to the scattering either. 

Concentrations shall be written as dimensionless 
mol/mol fractions, normalized to a pure NaAlH4 sam- 
ple. Accordingly, the initial concentrations are [A] = 1, 
[B] = [C] = 0. Since our scattering signal is almost ex- 
clusively due to bound hydrogen, we set 

aiW = \[B]fB, (6) 
a 5 (t) = [A]f A + \[C\f c , (7) 

where the / are Lamb-Mossbauer factors. 

On this base, we searched for rate equations that re- 
produce the observed time series ai (t) and a$(t). We first 
used the code generator kinpy [39j to try kinetic models 
based on stochiometric reaction equations. It turned out 
that the relatively sharp peak in a\(t) can only be re- 
produced if some autocatalysis is assumed. To improve 
the agreement with the experimental ag(t) and ai(t), we 
admitted concentration dependences with fractional ex- 
ponents. In this heuristic way, we found that the fol- 
lowing kinetic model is about the simplest one that is 
compatible with the measured as(t) and a\{t): 

d[A]/dt = -koa[A] -k 01 [A][B} 2 , 

d[B]/dt = -d[A]/dt - d[C]/d*, (8) 

d[C]/dt = fc 10 [£] 4/3 + fcii[£] 4/3 [C] 2/3 , 
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At 185°C, the following rate coefficients are found: 



fc 00 = 0.029 h~\ 

jfcoi = 0.55 h~\ 

fcio = 0.016 h _1 , 

fcn = 0.097 h" 1 . 



(9) 



The spontaneous rate coefficients kjo (j = 0, 1) are 
smaller by about an order of magnitude than the au- 
tocatalytic ones kji, and the second autocatalytic step 
is much slower than the first one, as was to be expected 
from the shape of <zi(t). 

For the fits at different temperatures in Fig. [9l the 
spontaneous rates kjo were fixed. Besides a trivial time 
shift, the only adjustable parameters were the two auto- 
catalytic rates kj\. Their fitted values show a plausible, 
weak temperature dependence. The fits are not perfect, 
but this does does not necessarily indicate a failure of 
our kinetic model; it could also point to the limits of ex- 
perimental reproducibility (sample preparation, thermal 
history). 

The exponents in ((5J) cannot be determined very accu- 
rately; it can only be said that they are compatible with 
multiples of 2/3, as expected for an autocatalytic reac- 
tion that takes place at a surface. We suggest that kjo 
stands for nucleation and kji for crystallite growth. 

Our model is perfectly compatible with the reac- 
tion pathways proposed in a density-functional study 
He}. Based on first-principles calculations, the migrating 
species were suggested to be A1H 3 and NaH vacancies. 
The lowest activation energy was found for AIH3 vacan- 
cies, leading to the following pathway for reaction (JTJ: 



nNaAlHi -> nNaAlH 



aih ;! 



Al 



2 H 2, 



(10) 



where the superscript denotes one vacancy, and n in- 
dicates an arbitrary amount of bulk material. The 
proposed diffusion mechanism also involves an AlHg~ 
ion. Anyway, the vacancy ultimately reaches a NaAlH4- 
NaaAlHg boundary, where it annihilates, releasing an 
excess Na + that aggregates with the growing NasAlHg 
phase. From another computational study we learn that 
this growth additionally requires a two-step transforma- 



tion A1HJ -> AlHg 



AlHfi \24 



Coming back to Ref. [40( , an alternate pathway starts 
with the unassisted release of Na + and H~ at the bound- 
ary of the growing Na 3 AlH 6 according to 

nNaAlH 4 -> (ra - l)NaAlHl NaIi + Na 3 AlH 6 . (11) 

The two ionic vacancies migrate together to the Al- 
NaAlH4 boundary, where hydrogen is released: 



nNaAlH^ aIi 



-> (n - l)NaAlH 4 + Al 



|H 2 . 



(12) 



Using positron annilation, it was confirmed experimen- 
tally that vacancies are formed in NaAlH4 during dehy- 
drogenation 41 1. 

In these pathways, boundaries play a crucial role: The 
NaAlH4-Na3AlH6 boundary as a sink for AIH3 vacancies 



and as a source for NaH vacancies, and the Al-NaAlH4 
boundary as a sink for NaH vacancies and as the location 
of H2 release. At least in undoped material, hydrogen 
release rates arc limited by processes at these boundaries. 
This is in accordance with the autocatalytic terms of our 
kinetic model. 



C. Effect of Doping 

In Fig. [TU1 the kinetic model ([5]) was applied to TiCl3- 
doped NaAlH4, cycled eight times under hydrogen and 
quenched in the hydrogenated state. These fits were at 
best a very rough first approximation, because constant 
rate coefficients were assumed although the sample tem- 
perature varied between 40 and 137°C (the excursion to 
167. . . 182°C was excluded from the fit). The model may 
even be entirely inappropriate, because the rate-limiting 
processes in doped alanate could be qualitatively differ- 
ent from those in pure material. According to [40], hy- 
drogen release in pure material is limited by surface pro- 
cesses, whereas in doped material rates are determined 
by vacancy diffusion. 

Under these reservation, the coefficients are about 



fcoo 
fcoi 
fcio 



0.05 h" 
0.8 h -] 
0.12 IT 



(13) 



fcn had no noticeable influence upon the time dependence 
of as(t) and ai(t). All three values in (fT3|) are higher than 
their counterparts in © although the latter were deter- 
mined at the considerably higher temperature of 170°C. 
Doping was expected to increase the coefficient fcoo, so 
that the first decomposition step starts at much lower 
temperatures than in the undoped material. However, 
the effect of doping on the second decomposition step 
seems to be even stronger: The coefficient kio is larger 
than in the undoped 185°C data by a full order of magni- 
tude. This correlates with the observation that the peak 
height of a\(t) is much lower in the doped than in the 
undoped material. 



VII. CONCLUSION 

The high flux, low background, and fine resolution of 
the new backscattering spectrometer SPHERES opens 
up new perspectives for time-resolved studies of weak and 
narrow quasielastic spectra. The present work reveals 
the potential of this method for investigations of reaction 
kinetics. 

The desorption of hydrogen from sodium alanate is an 
almost perfect test case, because on the one hand with 
its five different reactants it is complex enough to be 
interesting, while on the other hand the neutron spectra 
are not too complex, because two of the reactants do not 
noticeably contribute to the scattering, and only one of 
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them, the intermediate reaction product NasAlHg, shows 
quasielastic scattering. Once these quasielastic spectra 
are well characterized, the time-resolved measurements 
allow for an automatic data reduction that yields a small 
number of time-dependent amplitudes. These amplitudes 
can be related in a straightforward way [Eqs. ([6]), (0] to 
reactant concentrations. 

A time-resolved measurement of NaAlti4 decomposi- 
tion at 170. . . 185°C was fitted by a physically meaningful 
set of reaction rate equations. According to this model, 
both desorption steps (jXJ) , © are controlled by nucleation 
and growth. 

Next, it will be interesting to investigate more system- 
atically the temperature dependence of reaction rates. 
To initiate the decomposition of NaAlH4 , there seems to 
be only a small temperature range between about 165°C 
and the melting point. However, once the reaction has 
started, its autocatalytic sequel possibly remains active 
at considerably lower temperatures. 

While special care was taken to produce all samples in 
the same manner, it would be interesting to investigate 
how the kinetics depends on the initial texture and on 



the hydrogen loading/unloading prehistory. 

Most importantly, the basis is now available for sys- 
tematically investigating how the single reaction steps 
are modified by the presence of a dopant. As a very pre- 
liminary result, we hypothize that TiCl3 mainly acceler- 
ates the spontaneous decomposition of the intermediate 
reactant NaaAIHg. However, this still needs thorough 
verification. Unfortunately, only a very narrow quasielas- 
tic signal can be observed at the low reaction tempera- 
tures made possible by doping (Fig. [6]). Nevertheless, 
Fig. [TU1 suggests that a quantitative study of reaction ki- 
netics will be possible if more systematic time runs will 
be measured. 
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